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ABSTRACT 

We present the ultraviolet (UV) color-color relation of early-type galaxies (ETGs) in the nearby universe 
(0.05 < z < 0. 12) to investigate the properties of hot stellar populations responsible for the UV excess (UVX). 
The initial sample of ETGs is selected by the spectroscopic redshift and the morphology parameter from the 
Sloan Digital Sky Survey (SDSS) DR7, and then cross-matched with the Galaxy Evolution Explorer (GALEX) 
Far-UV (FUV) and Near-UV (NUV) GR6 data. The cross-matched ETG sample is further classified by their 
emission line characteristics in the optical spectra into quiescent, star-forming, and AGN categories. Contami- 
nations from early-type spiral galaxies, mergers, and morphologically disturbed galaxies are removed by visual 
inspection. By drawing the FUV -NUV (as a measure of UV spectral shape) vs. FUV - r (as a measure of 
UVX strength) diagram for the final sample of ^3700 quiescent ETGs, we find that the "old and dead" ETGs 
consist of a well-defined sequence in UV colors, the "UV red sequence", so that the stronger UVX galaxies 
should have a harder UV spectral shape systematically. However, the observed UV spectral slope is too steep 
to be reproduced by the canonical stellar population models in which the UV flux is mainly controlled by age 
or metallicity parameters. Moreover, 2 mag of color spreads both in FUV -NUV and FUV -r appear to be 
ubiquitous among any subsets in distance or luminosity. This implies that the UVX in ETGs could be driven 
by yet another parameter which might be even more influential than age or metallicity. 

Subject headings: galaxies: elliptical and lenticular, cD — galaxies: evolution — galaxies: stellar content — 
ultraviolet: galaxies 



1. INTRODUCTION 

The ultraviolet excess (UVX) refers to the phenomenon 
that the early-type galaxies (ETGs) exhibit a rising flux 
with a decreasing wavelength from 2500A to the Lyman 
limit. Theoretical and observational evidence shows that 
the phenomenon originates from the hot heli um-burning 
horizontal-branch (HB) stars and their progeny (lO'Connelll 
Il999t iBrown et alJl2000i) . One of the most important obser- 
vational constraints to unders tanding this phenomenon was 
given by iBurstein et alj d 1988b . They first explored the sys- 
tematics of the UV radiation from the 3 1 local ETGs and the 
bulge of M31, and found that the UVX strength correlates 
with the nuclear spectral line index Mg2 (and, weakly, with 
the central velocity dispersion and luminosity) in quiescent 
galaxies. With the advent of the Gala xy Evolution Explorer 
(GALEX) satellite dMartin et al.ll2005l) . a large area, deep UV 
imaging survey gathers high-quality data for an unprecedent- 
edly large sample of ETGs in the nearby universe. With the 
GALEX UV da ta cross-matched with the Sloan Digital Sky 
Survey (SPSS). [Rich et all (120051) investigated the systemat- 
ics of the ETGs at z < 0.2, but found no correlation of the 
UVX strength against the Mg2, D 4000 indices or stell ar veloc- 
ity dispersion. On the other hand, Donas et al] d2007l) found a 
tight correlation between the GALEX UV color and the Mg2 
index a mong the morphologically sorted ETG s in the RC3 
catalog Jde Vaucouleurs et all ll 991). Recently. iBureau et alJ 
(1201 ll) recovered the correlation of Burstein et al. from the 
GALEX UV imagings combined with the high-quality data of 
the SAURON sample, and also with the SDSS data for a more 
distant ETG sample. They argued that such correlations could 
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be seen among the quiescent ETGs only when the contamina- 
tions of star-forming galaxies and low-qualit y data were muc h 
more stringently removed than was done in lRich et alj (120051) . 

The UV color-color diagram is another powerful tool to in- 
vestigate the properties of hot stellar populations responsible 
for the UVX. Particularly, the capability of the GALEX to ob- 
serve the far-UV (FUV; 152 8A) and near-UV (NU V; 2271 A) 
bandpasses simultaneously (Morrissey et al. 2005) let us se- 
curely estimate the UV spect ral shape without in strumental 
biases. As demonstrated by Dor man et alJ (Il995h . the UV 
color-color diagram, FUV -NUV (as a measure of UV spec- 
tral shape) vs. FUV - r (or FUV -V; as a measure of UVX 
strength) can trace the UV spectral shape and amplitude. It is 
also usefu l to se gregate different stellar population systems. 
iRich et all (120051) employed this tool but f ound no system- 
atic variation among their ETGs at z < 0.2. iRee et all (2007) 
presented another example to show that the brightest clus- 
ter galaxies (BCGs) are homogenous populations following 
an old stellar sequence with considerable amounts of color 
spread. 

In this Letter, we expand significantly the sample size of 
previous studies from the latest data release of the GALEX 
and SDSS to construct the UV color-color relation of the 
ETGs in the nearby universe, and reexamine our current un- 
derstanding of UVX. 

2. GALEX - SDSS EARLY-TYPE GALAXY SAMPLE 

The initial sample of ETGs is selected from the SDSS DR7 
spectroscopic and photometric database in a rather simplistic 
way, with the following two criteria: 

spectroscopic redshift : 0.05 < z < 0. 12 
morphology parameter : fracDeV r > 0.99 

The search volume is limited as above in order to secure an 
unbiased sample from the SDSS spectroscopic survey. At 
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Figure 1. (Left) : u — r vs. M r color-magnitude diagram (top) and u — g vs. u — r color-color diagram (bottom) of the initial SDSS ETG catalog. (Right) : 
FUV -r vs. M,- color-magnitude diagram (top) and FUV —NUV vs. FUV — r color-color diagram (bottom) of the GALEX-SDSS ETG cross-matched catalog. 
The symbols are color-co ded according to the spectroscopic types (red: quiescent, green: AGN, cyan: composite, blue: star-forming) by the emission line 
diagnostics I Oh et al. 201 1). Filled circles with error bars indicate the median colors and magnitudes with the standard deviation of each classified subgroup. 
Some E+A galaxies are found in our data and their UV colors are indicated by asterisks. The median observational errors are denoted by the black crosshairs. 



z < 0.05, the SDSS spectroscopy is incomplete for bright 
galaxies. The SDSS target selection algorithm of the main 
galaxy sam ple allocates the fibers to the galaxies r = 17.77 
or brighter dStoughton eTaT1l2002t fStrauss et al.ll2002l) . This 
magnitude limit corresponds to M r w -21 at z = 0.12. Then, 
we measure the early-type morphology with the fracDeV 
parameter in r-band. The fracDeV is a fractional parame- 
ter that measures the amount of flux contribution by fitting a 
galaxy image to a de Vaucouleurs surface brightness profile 
(Ide Va ucouleurs! [T948h . By selecting fracDeV r ~ 1 galaxies, 
our sample secures purely spheroidal galaxies even though a 
few percent of the sample turn out to be disk galaxies with 
prominent bulges. Those contaminants are removed later 
from the final sample through visual inspection. The SDSS 
DR7 databas^3 returns an optical catalog of 94,796 galaxies 
with our querying parameters. 

The GALEX GR(fl provides the cross-matched table 
(xSDSSDR7) against SDSS DR7. We search this table for 
the UV sources (detected both in FUV and NUV images) 
matched against the input optical catalog and find a total 
of 34,388 UV sources cross-matched to a total of 23,493 
unique SDSS objectlDs. About half of the matched sources 
are observed in the GALEX All-sky Imaging Survey (AIS; 
texp ~ 110 sec), and the other half in the Medium Imaging 
Survey (MIS) of a single orbit exposure (t exp « 1500 sec) or in 
the Deep Imaging Survey (DIS) of multiple orbits. Due to the 
short exposure time, the GALEX AIS data are highly biased 

3 http://cas.sdss.org/dr7 
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toward star-forming galaxies, whereas the quiescent ETGs are 
well defined in deeper exposure images. The GALEX AIS 
data are inadequate for the demographic study of ETG popu- 
lations at these redshifts and hence will be ignored hereafter. 

For those ETGs detected in sufficiently deep exposure im- 
ages with the GALEX, we analyze their emission line charac- 
teristics of the optical spectra. The emi ssion line stren gths in 
the SDSS DR7 spectra are measured bv lOh et alj (120111) . and 
the sample galaxies are cl assified by the emiss i on line ratios, 
rOIIIl /H/3 and rNIIl/Ha dBaldwin et al.1 fT98~il [Kewiev et all 
I200U iKauffmann et al.1 120031) . lOh et alj d201 lb provide the 
amplitude-to-noise (A/N) ratio at each emission line to set 
the analysis threshold. As our main purpose is to secure the 
quiescent (i.e., no emission line) ETGs as stringently as possi- 
ble, we loosen the threshold to A /N = 1 and classify the sam- 
ple galaxies into quiescent (33%), star-forming (7%), com- 
posite (12%), and AGN (48%) categories. 

Fig. 1 shows the color-magnitude diagrams and color-color 
diagrams of the cross-matched ETG catalog. In the left pan- 
els, our sample of ETGs is expressed in the u, g, and r-band 
magnitudes and the colors sub-grouped according to the emis- 
sion line characteristics. The separation among the differ- 
ent spectroscopic subgroups (especially between the quies- 
cent and AGN-host galaxies) is much more effective in UV 
colors, as shown in the right panels. Most of the quiescent 
(red) galaxies are located in the "UV red sequence", whereas 
the blue envelop is occupied by star-forming (blue) galaxies. 
AGN-host (green) galaxies are spread out in-between in the 
UV color-color diagram. We adopt the modelMag (SDSS) 
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and MAG_AUTO (GALEX) as the total magnitudes of the 
galaxies in each band. Even though the source extraction 
apertures are not identically applied between the photome- 
try of the optical and UV images, the ap erture effect is sm all 
compared to the observa tional errors (see Ree et al. 2007). 

Some E+A galaxies (iGotol 12007) are found in our data 
and their UV colors are indicated by asterisks in Fig. 1. 
E+A galaxies are located over the wide color ranges (0 < 
FUV-NUV < 3, 2 < FUV-r < 8), but clearly separate 
from the "UV red sequence". They appear to be on an evo- 
lutionary sequence in which the UV flux fades away after the 
starburst, and yet FUV declines much faster than NUV (e.g., 
iKavirai et alJl2007at iChoi eFaIll2009h . 

3. UV COLOR-COLOR RELATION OF QUIESCENT EARLY-TYPE 
GALAXIES 

We examine the morphology of individual galaxies in 
the SDSS optical images for those 4025 spect roscopically 
quiescent ETGs. There are still some contaminants of early- 
type spiral galaxies, mergers, or morphologically disturbed 
galaxies in our sample. We find that about 9% of our sample 
should be rejected in the final sample of quiescent ETGs. This 
number also includes some possible blends in the GALEX 
UV images in the case that there are close neighbors of blue 
dwarf galaxies around a target ETG Fig. 2(a) shows the col- 
ors of such contaminants (denoted by crosses). Note that we 
eliminate them by the optical morphology only, being blind 
to their UV or optical colors. Most of those morphological 
contaminants enter the quiescent sample because their optical 
spectra are taken within the small fiber aperture (3" in diam- 
eter) toward the galactic centers, while they do contain young 
stars in the outer region judged from the strong extended 
UV emission. The remaining 3677 morphologically-cleaned 
quiescent ETGs consist of a well-defined "UV red sequence" 
as denoted by contours in Fig. 2(a), with the peak number den- 
sity at FUV-NUV ks 1.1 and FUV - r « 7.0. The UV col- 
ors of the final sample of quiescent ETGs at 0.05 < z < 0.12 
are very similar to the color distribution of the local ETGs 
(Fig. 2b). The GALEX RC3 E (filled circles) and SO (trian- 
gles) galaxy data are from Donas et al. (2007), and their V 
magnitudes are converted to r-band with V -r = 0.29 assum ed 
from the 12 Gyr elliptical galaxy model in Ree et al.1 (120071). 

Base d on the semi-empirical criterion of Schawinsk i et al.l 
d2007l) . 31% of the quiescent ETGs in our sample are classi- 
fied as the recent star formation (RSF) galaxie s (NUV — r< 
5.4), w hich is consistent with the f indings of iKaviraj et al.1 
(2007b|) and lSchawinski et al.l (120071) . Note, however, that the 
observed UV colors of quiescent ETGs follow the sequence 
of a much steeper slope rather than the semi-empirical model 
line (i.e., a harder UV spectral shape at a given UVX am- 
plitude). Using a sigma-clipping linear fit algorithm to the 
observed colors iteratively, we derive the color-color rela- 
tion, (FUV-NUV) = 1.36 (FUV-r)-&35, for the quies- 
cent ETGs, which indeed has a steeper slope than that of the 
semi-empirical model line, i.e., NUV — r = 5.4. In Fig. 2(c), 
we empirically define the "UV red sequence" (red dots) from 
the linear fit to the observed colors (solid line) with 1 .5 mag 
color width (which corresponds to 3 times the mean absolute 
deviation of the observed data from the fit line) in FUV-r 
(dotted lines). The UV colors of the BCGs in Abell clusters 
(filled circles) fall well into this sequence. There are still a 
significant number (20%, blue dots) of outliers from the UV 
red sequence with this empirical definition, and their addi- 
tional UV flux should be attributed to the presence of small 
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Figure 2. FUV-NUV vs. FUV — r color-color diagram of the quiescent 
ETG sample, (a) : Among 4025 spectroscopically quiescent ETGs, about 9% 
of the sample (denoted by crosses) turn out to be early-type spiral galaxies, 
mergers, or morphologically disturbed galaxies in the visual inspection of 
SDSS images. The remaining 3677 morphologically cleaned ETGs consist 
of a well-defined "UV red sequence" (denoted by contours), (b) : The UV 
colors of the final sample of quiescent ETGs at 0.05 < z < 0.12 are com- 
pared to the color distribution of the RC3 E (filled circles) and SO (triangles) 
galaxies, (c) : The "UV red seqence" is empirically defined by the linear fit 
(solid line) with 1.5 mag c olor width in FU V — r (dotted lines). The colors of 
Abell BCGs (filled circles. iRee et alj 2"0071 fall well into t he sequence. The 
semi-e mpirical criterion of RSF galaxies (NUV —r = 5.4) in Schawinski et al. 
12007) is indicated by the dashed line for comparison. 

fractions (~1% in mass) of young stars (lYi et alJl2005l) . The 
number fraction of such outliers is negligible (< 5%) in the 
brightest (M r < -22) subgroup, but gradually increases to- 
ward the less bright galaxies (~ 30% forM,. » — 21 subgroup). 
We find no substantial difference between the lower redshift 
(0.05 < z < 0.086) and the higher redshift (0.086 < z < 0.12) 
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subgroups. 

All the colors of galaxies presented here are derived 
from the apparent total magnitudes in each filter with no 
other corrections applied, except for th e Galactic extinction 
(ISchlegel et al.lfl998h IWvder et alJl200l . The effects of red- 
shift (0.05 < z < 0.12) and evolution (At « 1 Gyr) are small 
compared to the intrinsic color spreads among the sample 
galaxies and should have little impact on the general mor- 
phology of the color-color diagram. However, better statistics 
(e.g., the number fraction of RSF galaxies) require the sam- 
ple galaxies to be at equivalent cosmological distances. More 
detailed demographic analysis is in preparation and will be 
presented in a future paper. 

4. MODEL COMPARISONS 



We employ the x 2 minimization technique of lYi et alJ 
(1201 ll) to find the best fit model spectra for the observed multi- 
band photometry. Most of the observed data in Fig. 2 are well 
fitted by the two-component models in which the UV spectral 
shape is determined by the combination of the baseline old 
(12 Gyr) stellar spectra and the varying amount (0.01 - 100%) 
of a young (0.01 - 10 Gyr) stellar component. However, the 
fitting accuracy is poor for the sample of quiescent galaxies 
in the lower right corner (FUV-NUV < 0.9, FUV-r > 6.0) 
of Fig. 2. The galaxy in this particular color domain has a 
very hard UV spectral shape (comparable to that of NGC 1 399 
which is one of the strongest UVX galaxies in the local uni- 
verse), but its UVX amplitude is relatively small. 

The observed UV colors are compared with other popu- 
lation synthesis models in Fig. 3. We first adopt the com- 
po site stellar popul ation (CSP) model (t = 12 Gyr, Z = Z©) 
of iMarastoril d2005l) for the baseline model of bulk popula- 
tions which does not contain blue HB stars. Then, in or- 
der to simulate the UV color variations due to different types 
of UV sources, we add the simple stellar populations (SSPs) 
of Z = 0.1 Z Q for the young components (t = 0.1, 0.5, and 
1.0 Gyr) and for the old compone nt (t = 12 Gyr, 100% blue 
HB) from IConrov & Gunnl d2010l) models. The recent mod- 
els of lChung et alJ d201 ll) are also considered for the helium- 
enhanced cases (Y = 0.33, 0.38, and 0.43). The flux weight of 
SSPs is assumed to be 10% of the baseline CSP model in r- 
band. This weight translates roughly to 0.3% (0.1 Gyr), 1.5% 
(1.0 Gyr), and 10% (12 Gyr) in stellar mass fraction. The 
top panel shows the UV-optical model spectral energy distri- 
butions (SEDs) for the various possible UV sources, and the 
UV colors of each toy model are derived from the CSP+SSP 
model SEDs redshifted to z = 0.086 (median redshift of our 
sample) as shown in the bottom panel. 

The shaded region in the bottom pan el of Fig. 3 denotes the 
canonical models in iRee et all d2007l) where the mean tem- 
perature of HB stars smoothly varies accor ding to metallic - 
ity and age dLee et al.l[l994t lYi et all 1 19991) . It is clear that 
the observed color spreads are much larger than the canon- 
ical model predictions spanning large ranges of age (11 - 
15 Gyr) and mean metallicity (0.5, 1.0, and 2.0 Z ), and the 
difference cannot be attributed to observational errors or red- 
shift effects. Moreover, the observed slope of the color-color 
relation is much steeper than the model predictions, imply- 
ing that the UVX sources should have a harder UV spec- 
tral shape (i.e., hotter temperature) at a given amplitude sys- 
te matically. An ad hoc assumption in the S EP combination 
oflMarastonl d2005l) CSP (no blue HB) and IConrov & Gunnl 
d2010l) SSP (100% blue HB) models could produce a some- 
what harder UV spectral shape at a given UVX amplitude 
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Figure 3. (Top) : UV-optical model SEDs. As the Marastorj <2005T) CSP 
model (12 Gyr, solar metallicity) does not contain blue HB stars, the SSP 
models (Z = 0.1 Zq) for young components (t = 0. 1, 0.5, and 1.0 Gyr) an d 
for old component (t = 12 Gyr, 100% blue HB) from lConrov & Gunnl (20101) 
and for helium-enhanced cases (Y = 0.33, 0.38, and 0.43) from Chung et al. 
1 201 1 ) are additionally adopted to reproduce the UV spectra. The flux weight 
of SSP models is assumed to be 10% of CSP in r-band. (Bottom) : UV colors 
derived from the model SEDs redshifted to the median redshift of the sample. 
The combinations of CS P+SSP SEDs are as labeled. The canonical stellar 
population model grids in Ree et al. ( 2007) are denoted by the shaded region 
(t = 1 1 - 15 Gyr, < Z > = 0.01, 0.02, 0.04). 



(e.g., FUV-r rj 7) than the models in lRee et all d2007l) . and 
yet it is still hard to explain the UV spectral shape of the 

strongest UVX galaxies at FUV - r 6. 

We test the recent model prediction of IChung et alJ (1201 ll) 
to determine whether the observed color spreads and the steep 
slope of the color-color relation might be reproduced with 
the helium-enhanced stellar populations. Since the helium- 
rich stars have lower masses at a given age, the mean tem- 
perature of helium-rich HB stars is much higher than the 
normal-helium populations. This effect is demonstrated well 
in Fig. 3. The combinations of the baseline CSP model 
and the helium-enhanced (Y = 0.33 and 0.38) SSP models 
would probably explain the UV spectra of the strongest UVX 
(FUV-r w 6) galaxies under the ad hoc assumption that all 
the UV flux originates from the helium-enhanced HB stars 
only. An even higher helium abundance (e.g., Y = 0.43) seems 
unlikely to contribute in this color-color diagram as the HB 
stars of such high helium abundance are too hot and faint in 
the GALEX FUV bandpass. The UV spectral shapes of the 
helium-enhanced stars are clearly disentangled from those of 
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the young stars which have a much stronger UV flux both in 
FUV and NUV at the same optical luminosity. 

5. CONCLUSION 

In summary, we made two main findings with the latest 
GALEX UV data of the ETGs in the nearby universe. 1) The 
"old and dead" ETGs consist of a well-defined sequence in 
UV colors in which the stronger UVX galaxies should have 
a harder UV spectral shape systematically. The UV color 
distributions of ETGs at 0.05 < z < 0.12 are very similar to 
those observed for the local RC3 galaxies. The comparison 
with the population synthesis models indicates that the hot 
HB stars are responsible for the observed "UV red sequence", 
while the contamination of recently-formed young stars can 
be effectively discriminated in the UV color-color diagram. 
2) The observed slope of the UV color-color relation of qui- 
escent ETGs is too steep to be reproduced by the canonical 
stellar population models in which the UV flux is mainly con- 
trolled by age or metallicity parameters. Moreover, 2 mag 
of color spreads both in FUV -NUV and FUV — r appear to 
be ubiquitous among any subsets in distance or luminosity, 
which can hardly be explained by the canonical models with 
even large spreads in age (> ±2 Gyrs) and mean metallicity 
(> ±0.3 dex) assumed. This implies that the UVX in ETGs 
could be driven by yet another parameter which might be even 
more influential than age or metallicity. 

The recent debate on the helium-enhanced populations has 
an interesting implication for this issue. Observations sug- 
gest that there are some helium-enhan ced subpopulations 
in th e extend e d HB globular c lusters (D'Antona & Calol 
20041 iNorrisI l200l [Lee et al.1 |200H JPiotto et al.l l2005t 
Recio-Blanco et al.ll2.006t iLee et al.ll2007l) a nd in the massive 
globular clusters in M87 (ISohn et al.ll2()06). If present, such 
subpopulations would steepen the UV spectral slope system- 
atically by increasing the mean temperature of helium burn- 
ing stars, as demonstrated in Fig. 3. They would also en- 
hance the UV color spreads among the ETG stellar popu- 
lations with even small variations in the numb er fractions 
or in the amount of helium enhancement ("e.g.. iKavirai et al] 
I2007ct IChung et all 120111) . Whether the postulated helium 
enhancement in old stellar systems would occur also in galac- 
tic scales or in the intracluster medium of galaxy clusters 
dPeng & Nagail |2009|) is in que s tion. Obs ervational tests 
by lLoubser & S anchez-Blazquea d201 ll) and lYi et al.l (1201 lh 
do not supp ort the helium sedimentation hypothesis, while 
iCarter et al.ld20Tl argue that the helium abundance is a plau- 
sible candidate to explain the UV radial gradient. By the com- 
parison between the UV color distribution and some toy mod- 
els presented in this Letter, it is suggested that the observed 
steep slope of the UV color-color relation of ETGs would 
be better explained with the helium-enhanced subpopulations. 
The excessive spreads in the UV colors observed among the 
ETGs, besides the variations expected by the age/metallicity 
spreads, would be attributed to small variations in helium 



abundance or to the low-level star formation that is hard to 
be identified with the current survey depths and resolutions of 
the GALEX and SDSS data. 
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